Toroidal Perturbations of Friedmann-Robertson- Walker Universes 
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Explicit expressions are found for the axisymmetric metric perturbations of the closed, flat and 
open FRW universes caused by toroidal motions of the cosmic fluid. The perturbations are decom- 
posed in vector spherical harmonics on 2-spheres, but the radial dependence is left general. Solutions 
for general odd-parity /-pole perturbations are given for either angular velocities or angular momenta 
prescribed. In particular, in case of closed universes the solutions require a special treatment of the 
Legendre equation. 

PACS numbers: 04.20-q, 98.80.Jk 



I. INTRODUCTION 

There exists a class of perturbations in the Newtonian 
astrophysics which represent differential rotations of a 
spherical star in hydrostatic equilibrium. The displace- 
ment fields are time-independent and have no radial com- 
ponents Q: 



4wGp£ r = 0, 4ivGp& 
4irGp£ v -- 



Timjr) dYi m 
r sin 9 dip 
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r 80 
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where Ti m is any function satisfying the condition that 
£e and £ v are finite. These solutions define the toroidal 
perturbations. When 1 = 1, the displacements can be 
written as 



£ = 5Sl(r) x r, 



\sn\ 
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which represent rigid rotations of the internal spherical 
surfaces around the center. If <5f2 docs not depend on 
r, the whole star is rotated as a rigid body. To get a 
uniform slow rotation, rather than a static displacement, 
we assume equation (ll.ljl to be true with £ replaced by 
the Eulerian velocity field; then displacements are linear 
in time. Such perturbations are analogues of the cosmo- 
logical toroidal perturbations investigated in this paper. 

The toroidal perturbations (|l.lf> are transverse and di- 
vergence free. They are often called the "trivial modes" . 
The trivial modes become important when one considers 
the second order perturbations (~ [<5ri] 2 ) 2]. They de- 
velop into the so-called r-modes in non-radial oscillations 
of a rotating star (see e.g. 3]). They also give rise to non- 
zero eigenvalues when an external (e.g. magnetic) field 
is added. 

Toroidal oscillations are important in geophysics 0. 
Non-rotating, spherically symmetric Earth models ad- 
mit "trivial toroidal modes" which are associated with 
vanishing eigenfrequencies; they do not alter the elastic- 
gravitational potential of the Earth. As in astrophysics, 



their counterparts play a significant role on a rotating 
Earth. In the solid core the toroidal modes are nontrivial 
even in spherical models: the presence of an anisotropic 
stress tensor produces a traction. 

The velocity pattern of the axisymmetric I = 2 toroidal 
(quadrupole) mode is shown in Fig. 1. The individual 
rings of fluid rotate independently, preserving their local 
angular momentum. For non-trivial modes in the solid 
core the displacements will reverse after a half of the 
oscillation cycle. 




FIG. 1: Velocity field for axisymmetric quadrupole (I — 2) 
toroidal perturbations. 

In general relativity even "trivial modes" in spherical 
systems become "non-trivial" since they produce drag- 
ging of inertial frames and, hence, influence essentially 
gravitational field. 

In the following we investigate the toroidal perturba- 
tions of Friedmann-Robertson- Walker (FRW) universes 
of all three types (k = — 1,0, +1). We confine ourselves 
to the axisymmetric case. Since the backgrounds ad- 
mit homogeneous and isotropic foliations, non-symmetric 
perturbations can be found from the axisymmetric ones. 

Indeed, in any spherically symmetric background sys- 
tem, one can restric oneself to axisymmetric modes of 
perturbations since non-axisymmetric modes with an 
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e %mip dependence (m = —I, . . . , +1) can be derived from 
the axisymmetric perturbations (m = 0) by suitable ro- 
tation of the axes. If, in the original coordinates, the new 
polar axis will be pointing in a direction (9' ,ip'), then an 
axisymmetric mode with respect to the new axes can be 
decomposed in non-axisymmctric modes in the original 
axes by using the addition theorem for spherical har- 
monics (see 0, P- 138-139 for more detailes). Each 
the e ml(p component will separately solve the perturba- 
tion equation with the original radial functions. Since 
we are studying vector perturbations, the vector compo- 
nents will also transform under a rotation. (So, in (|1.1|) 
with an axisymmetric mode only £ v ^ 0, but a general 
rotation will induce £g ^ as well.) 

The equations to be solved are derived and discussed in 
Section 2 of |5j], hereafter Paper I. Here we just summa- 
rize few basic relations. We start out from the perturbed 
FRW metric in the form 



ds 2 = (g^ + h^dx^dx" 

= dt 2 - a 2 (t)f i jdx i dx j + h^dx^dx" , 



(1.3) 



where the background metric ~g is used to move in- 
dices; the time-independent part of the spatial back- 
ground metric [i,j,k = 1,2,3] is used to define the 
3-covariant derivative Vfc, V fc = / fci V;. We choose spa- 
tial harmonic gauge conditions on t = constant slices. 
The slices themselves are chosen so that the perturba- 
tion of the external curvature vanishes (uniform Hubble 
expansion gauge) - see Paper I for details. In addition, 
since we are interested in the toroidal perturbations we 
also assume the perturbation of the fluid velocity to be 
transverse, \?kV k — 0, and the same for the metric vector 
perturbation hok, V fe /iofe = 0. 

In the axisymmetric case we get the perturbed con- 
straint equation 



W 2 h 



Oip 



2kh Qv> = 2a 2 K6T° 



(1.4) 



in which V 2 = / fci VfcV; and the perturbed Bianchi iden- 
tities (the dot denoting d/dt) 



(p + p)(a 2 r 2 sin 2 6V V - h 0ip )] = 0, (1.5) 



equivalently 



[a 3 5T°] 



0. 



(1.6) 



These just express the conservation of angular momen- 
tum of each element of the axially symmetric ring of fluid, 
analogously to the "trivial modes" in astrophysics men- 
tioned above. To have solutions in more closed forms, we 
decompose quantities only in vector spherical harmonics 
in variables 9, ip, but leave the radial dependence general. 
Hence, we write 



hotp = a 2 r 2 ^w;(t, r)sm9Yio,e, (1.7) 
i=i 

oo 

V v = -a 2 r 2 5^n,(t,r)flin0Y M , fl , (1.8) 



and 



5T° = a 2 (p + p)r 2 J2(ui -Cli)sm0 Y m>8 
1=1 

oo 

= 5}JZ£(t,r)]isiix0y,o ie . (1.9) 



i=i 



From equation l|1.4f> then follows the radial equation for 
each I in the form 

d 



r z or 



Vl - kr 2 — (r 2 uji) 
or 



l l±^-u; l -4kcj l = X 2 {n l -cj l ), (1.10) 

where A 2 = 2na 2 (p + p) = 4(fc — a 2 H), H — a/a is the 
Hubble constant; the last relation follows from the back- 
ground equations for arbitrary p, p, k and cosmological 
constant A. 

For k — equation (|1.10(1 can be written with the 
angular momentum density (6T®)i as a source in the form 

1 d ( 4 0wA 1(1 + 1) -2 a 



(1.11) 



for the fluid angular velocity as a source, 



1__5 

r 4 dr 



A2 + fcl) 



-A 2 a 



(1.12) 



here A 2 = -4a 2 H = 2na 2 (p + p). 

For k = ±1 it is advantageous to write r 2 = k(l — p 2 ), 
or p = Vl — kr 2 , to obtain 
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[k(l- p 2 )] 3 / 2 dp 



(t,M 2m5/2 5w M 1(1 + 1) -2 



dp 

2 k 



If we set 



k(l - p 2 ) v 

^ = [fc(i- M 2 r 3/4 



k(i-p 2 ) 

(1-13) 



(1.14) 



equation l|1.13(l becomes the Legendre equation for u>i 
with (6T®)i as the source: 



d_ 

dp 



kil-p 2 )^' 



duji 
dp 



fc 3 ( 3 d + W 

2 ( 2 +iJ k(l-p 2 ) 
2k 



[k(l - p 2 



W S 5T > (1-15) 



With the fluid angular velocity Qi as the source, the equa- 
tion reads 



d_ 

dp 



k{1 -^-p 



kv(v + 1) - 



(1 + 



1)2 1 

2> 



1 = 1 



k(l-p 2 ) 

K^-XMlkil-p 2 )} 3 ^, 



(1.16) 
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where 



i > + - ] = 4-2kna 2 {p+p) = 4-fcA 2 = Aka z H. (1.17) 



In the following we shall first solve equations l|l.llfl and 
(|1.12|) for flat universes (Section 2). The interesting but 
more complicated case is that of closed universes, i.e., 
solutions of equations (|1.15|) and l|1.16fl with k = +1. 

We shall discover that Legendre functions Q„_|/ 2 which 
solve equation (|1.15|) and are regular at least at one pole 
on the 3-sphere, in fact vanish identically for I > 2! How- 
ever, appropriate solutions can be found by considering 
the derivatives of the Legendre functions with respect to 
their degree v. We did not find this fact in the litera- 
ture on special functions. We analyze solutions of both 
equation l|1.15fl with the angular momentum given and 
(|1.16|) with the fluid angular velocity given for closed uni- 
verses in Section 3. Finally, in Section 4, we study them 
in detail also for open universes. Roughly speaking, the 
dragging effects are not damped when angular momenta 
are given as their sources, whereas they are, if fluid's an- 
gular velocity is considered as a source. The physical 
explanation of this apparent paradox is the same as that 
given in Paper I for I — 1 perturbations. 

Before we turn to solving the equations, we wish to 
emphasize that all our solutions are true at any given 
instant of time. However, the angular momentum and 
angular velocity distributions are related to those at ear- 
lier times by equations of motion (i.e., by the contracted 
Bianchi identities). In our case of toroidal perturbations 
they lead back to local conservation of angular momen- 
tum density, which demonstrates the advantage of taking 
angular momentum as the source. The time evolution of 
the dragging is thus independent of Z. It is analyzed in 
Section 5 of Paper I. 



II. AXISYMMETRIC TOROIDAL 
PERTURBATIONS FOR k = 

A. The angular momentum as a source of ui 

The homogeneous equation corresponding to (|1.11|) 
with {5T®)i given is 



4 , 1(1 + 1) -2 

U l + -UJ l ^ w i = 



(2.1) 



This can be easily solved by writing lui = r a ZU;, and 
choosing a so that the coefficient of uJi vanishes. It im- 
plies a = I — 1 or a = — (1+2) and the first-order equation 
for Two simple integrations lead to the solutions 



(|1.11|) which decays at infinity and is well-behaved at 
the origin is given by 



-2k 
21 + 1 



+ r 



(rT l (ST%dr' 



(2.3) 



In particular for the dipole (1 = 1) perturbations 



16tt 



' / V) 2 (£T°) 1 dr'+ r^ST^dr' 



If we write 



ST ( > = (£Z£) 1 sm0(-v / 3747rsin0), 



(2.4) 



(2.5) 



and integrate over 9, we find that l|2.4|l coincides with 
the second equation in (1.3.4); herafter, we denote in this 
way equations in Paper I. 

Equation 1)2. 3fl demonstrates how for general Z-pole 
distributions, "local" angular momenta, determined by 
(5T°)i, at all distances contribute instantaneously to the 
dragging of inertial frames without any exponential cut- 
off or influence of any horizon - see the detailed discussion 
of this point for I = 1 perturbations in Paper I. Notice 
that angular momenta (ST®)i for / > 1 do not, however, 
contribute to the total angular momentum in a spherical 
layer (r x ,r 2 ) given by [cf. (1.2.24) and (1.2.25)] 



J(ri,r 2 ) 



-2tt 



dr / d9 a 3 r 2 sin96T°. (2.6) 
/ i Jo 



Substituting for ST® from equation (|1.9fl . integrating by 
parts over 9 and realizing that 

/ Yio cos 6 sin 9d6 = ^4%/3 6 a (2.7) 
Jo 

due to the orthogonality of spherical harmonics, expres- 
sion l|2.6[) reduces to 



J(r u r 2 ) = 4^/3 a 3 drr 2 (5T°) l=1 

Jri 



(2.8) 



= A^/3a 5 (p+p) / drr 4 {uo l=1 -Q l=1 ). 
J r± 

Such result is understandable on intuitive grounds: e.g 
for I = 2 odd-parity perturbations rotations at antipodal 
9's are equal and opposite - the "northern" hemisphere 
rotates in the opposite direction to the "southern" hemi- 
sphere (see Fig. 1). 



B. The fluid angular velocity as a source 



Wl W =V-\ u\ IJ) =Bir-« +2 \ (2.2) 

A l ,B l constants. The Wronskian is W ; = AiB t (2l + 
l)r~ 4 . The solution of the inhomogeneous equation 



The equation to be solved is l|1.12|l . We introduce z 
Xr and assume A > 0, so z > 0. Putting then 



LJl = Z 3/2 'UJl, 



(2.9) 



4 



the homogeneous part of l|1.12fl turns into the equation For the source fli at z' ^> 1 we expand K l+ i [see l)2.12[l ] 



d 2 ^ 
dz 2 



1 duii 
z dz 



1 



1)21 



Ui = 0. 



(2.10) 



This is the modified Bessel equation, the general solution 
being given by 

oJ t = C l I l+i (z)+D l K l+t (z), (2.11) 



where Ii+i and -Kj+i are modified Bessel functions. For 
an integer / these functions are given explicitly by finite 
sums (see e.g. 8.467, 8.468 in 0, for still shorter expres- 
sions in terms of (d/dz) 1 , see 0, 0) which imply the 
following asymptotic expansions at z — > oo, 



K l+h {z) 




(2.12) 



so, independent of I in the leading terms at infinity. At 



the origin z 



L i+ 



0, 

2 1 
^(2/ + l)!l" 

7r(2/)! 1 

2~2~n\ s+i 



2(2/ + 3) 



z 2 + 



(2.13) 



The general solution of the original homogeneous equa- 



tion is thus given by loi 



u\ IIS> , where 



,U) 



A l I lH (z) = A l z-^ 2 I l+h {z), (2.14) 
B t K l+ i(z) = B lZ - 3 / 2 K l+ i(z), (2.15) 



be- 



(7;_)_i and are defined by these relations). 

haves well at z — ► 0, w[ at 2 — > 00. It is easy to see that 
for / = 1 one obtains the asymptotic expansions given in 
Paper I. The Wronskian reads (for the Wronskian of L , 1 

and K l+ i see [lfij . p. 68) 



(2.16) 



Using again the method of variation of parameters to 
find the solution of the inhomogeneous equation (|1.12|) 
which is well behaved at z — > and vanishes at z —* 00, 
we arrive at 



+ I /+ .(z) / {z'fK^iz')^')^', (2.17) 

and K l+ i are defined by Ij2.14jl and [|2.15jl . 
The solution (|2.17|l for general / can be explored in a 
way similar to that for I = 1, as was done in Paper I. So, 
to find ui for small z with sources f2; not close, we neglect 
the first term in H2.17J1 and expand in the second 

term: 



Vl(z) 



1 



tt(2/ + 1)!! 



1 



1 



2(2/ + 3) 



z 2 + 



L w 



K^iiz'^Wdz'. 



(2.18) 



1 



(2/ + 1)!! 



2(2/ + 3) 



z 2 + 



x / «) 2 e-* n,(y)d2'. 



(2.19) 



If, in addition, the source is localized in a small interval 
r (l±A), A< 1/A, we find 



1 



r(Ar) 



1-1 



(2/ + 1)!! 

x (Ar ) 3 e~ Aro H 2A 



1 



1 



2(2/ + 3) 



(Ar) 2 



(2.20) 



This generalizes equation (1.4.7) for any /. Thus, we 
recover the exponential decay which is now seen to be 
independent of /, but we also discover the power-law de- 
cline, given by (Ar) i_1 /(2/ + l)!! (Ar <C 1) which increases 
with increasing /. When f2; is concentrated near zq, in 
zq ± A A, then close to zq 1 and with Hi — £1/ we find 
wi(zo) = AAf^, as in Paper I. 



III. AXISYMMETRIC TOROIDAL 
PERTURBATIONS FOR k = 1 



We wish to solve the inhomogeneous equation (|1.13|) 
for uoi which, after the substitution i|1.14fl . turns into the 
inhomogeneous equation corresponding either to (|1.15|) 
or fTTfy 

So, to find solutions we have first to explore appro- 
priate solutions of the homogeneous Legendre equation 



d_ 

9/i 



cJl = (3.1) 



for / = 1,2,3,- •• and either v = |, or a general v 
given by H1.17|) . Since for the closed universe the variable 
/1 = cosx is real and /1 S (— 1, +1), we are interested in 
l+- i+- 

Legendre functions P v 2 (/z), Q v 2 (fx) on the cut. (We 
emphasize here that the order fj, of the Legendre func- 
tion, e.g. in has nothing to do with the inde pen dent 
variable fi.) Many useful formulas are given in 101. In 
particular, we start from the general expression for 
in terms of a particular combination of two hypergeo- 
metric functions (given as the 4th expression on p. 168 
in [T(il|'). change 6 — > x and put /i = / + i. Employing 
then the well-known properties of hypergeometric func- 
tions for special values of parameters, we find, after some 

algebra, that functions Q v 2 (cos x) can be written as fi- 
nite sums of the form 
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QL +2 (cosx) = 



i r(| + ^ + z) ' (-/)*(/ + i) fc / -i 



sin 1 / 2 x 



r(f 



E 



fc=0 



v 2 sin x , 



,(3.2) 



where the standard notation (a)k = cx(a + 1) • ■ • (a + 
fe — 1) = r(a + fe)/r(a) is used. In [lg, p. 168, suit- 
able expression in terms of hypergeometric functions is 
given also for Legendre functions P£ on the cut which 
one may wish to use as independent solutions. However, 
after putting /j, = I + i and making arrangements similar 

i+- i+- 
to those above for Q v 2 , it turns out that P v 2 diverge 

at x — 0, 7r. [For I = 0, 1 this was shown in Paper I.] 

j i i 

Therefore, only functions Q v 2 given by l|3.2|) will be ap- 
propriate objects in the following. 

It is easy to check that (|3.2H for I = implies equation 
(1.4.9). A small calculation shows that l|3.2[l gives 



i 

2 



sin 3/2 X 



(3.3) 



| — (n+ 1) sin x cos(nx) + sin [(n + l)x] j, 



which, after simple rearrangements, coincides with the 
first expression in equation (1.4.12) for 1=1. As in 



Paper I, we have put v = n — | in 1|3.3|1 and will do 
so hereafter. 

With (8T®)i given and k = +1, the homogeneous equa- 
tion corresponding to 11.15fl becomes (|3.1|) with v = |, 
n = 2. We thus will consider solutions Q 2 1 , n = 2, de- 

__ ™~ 2 

termined by (|3.2|l . In Paper I these solutions have been 
used for I — 1 and general n to determine u>i—i when Q/ = i 
is considered as a source. Since we also wish to consider 
solutions u>i for f2/ as a source with any Z, we do not put 
n = 2 yet, rather give first explicitly few functions 



Qn = Ql + \ 



(3.4) 



for general n and I = 1, 2, • • • , using the expression (|3.2|l . 
As with Z = 1, we shall also consider n purely imaginary, 
n = iiV. Calculations become lengthy with Z increasing. 
We have used MATHEMATICA for checking, simplifying 
and deriving some of the following formulas. 
Starting first from (|3.2|) for Z = 1, we get 



7T ~ 1 

Qn = ( 7T ) 2 . c, /o [cos X sin(nx) - » sin x cos(nx)] , (3.5) 
v2/ sm^^x 

C<w = * ) 2 . 3/2 [cos x sinh( A X ) - A sin X cosh( A X )] . (3.6) 
v 2/ sin^x 

For Z = 2 we obtain 

Qn = (|) 2 — ^-|-^sin(2x)cos(nx) + [2-n 2 sin 2 x + cos(2x)]sin(nx)|, (3.7) 

QiJV = *(f) 2 — ^|-^Asin(2x)cosh(Ax) + [2 + A 2 sin 2 X + cos(2x)]sinh(Ax)|. (3.8) 

With Z = 3 in (|3.2I) lengthier calculations lead to 



Ql = (|) 2 gin7/2 ~ ^ {" - 11) cos(2 X ) + n 2 ~ 19] cos(rcx) sin X 

+6 [(n 2 + l)cos(2x) - n 2 + 4] cosxsin(nx)}, (3.9) 
Q!n = *(f) 2 — ^i{A[(A 2 -ll)cos(2x)-A 2 -19]cosh(Ax)sinx 

+6 [(-A 2 + l)cos(2x) + A 2 +4] cosxsinh(Ax)}. (3.10) 



6 



Finally, explicit expressions implied by 13.2fl for I = 4 are 

Qt = (~Y . I, ^20n[(n 2 + 5)cos(2 X )-n 2 + 16]cos(nx)sin(2x) 



Q, 



JV 



(3 



/ sin 9 / 2 X 8 v 

3(n 4 - 25n 2 + 144) + 4(n 4 - 10n 2 - 96) cos(2x) - (n 4 + 35n 2 + 24) cos(4x) sin(nx) j , 
1 1 ' 



(3.11) 



20A [(—TV 2 + 5) cos(2x) + N 2 + 16] cosh(iVx) sin(2x) 



dn 9 / 2 x 8 

-3(A 4 + 25iV 2 + 144) + 4(7V 4 + 10A 2 - 96) cos(2x) - (N 4 ~~ 35iV 2 + 24) cos(4x)j sinh(TVx) \ ■ (3.12) 



The behavior of these solutions at the poles x = 0, 7T, as 
well as the second independent solutions to the homoge- 
neous Legendre equation i|3.1f) . will be discussed below. 

Now we have to distinguish the case when (ST®)i or £1; 
is considered as a source. 



A. The angular momentum as a source for cj 

In order to construct solutions to (|1.15|) for k = +1, 
we need solutions of the homogeneous Legendre equation 

(THl) for I = 1,2,3,... andi/= §, n = 2, i.e. Q?, = Q^. 

_ 2 
For Z = 1 we have the solution (|3.3I) with n = 2: 

Q2 = Q| =2(|) 5 sinix- (3-13) 

This is well-behaved at both x = and x = tt, whereas 
the second (independent) solution 



PI 



(I 



2 cos x 
sin^ x 



(l + 2sin 2 x) 



(3.14) 



diverges at these poles. The solution Q\ implies [cf. 
equation 11.1411 ] lui = \ = luq = constant, independent of x, 
as the solution of the original equation associated with 
(|1.13|) . Such wo may be time-dependent and may always 
enter u> in the case of closed universes, in which an "abso- 
lute rotation" does not occur. It cannot be removed (see 
equation (1.3.8), and text below). The other homoge- 

3/2 

neous solution for wi , corresponding to P 3 L , diverges at 

the poles. For I — 1 the inhomogeneous equation (|1.15|) 
can then be integrated directly. The integration yields 



du>x 
dx 



2k 



sin x Jo 



r S in 2 x'(«; )idx'. 
Jo 



(3.15) 



Starting from the definition of the angular momentum 
J(xi, X 2) m (1-2.25), with r = sinx, and integrating over 
8 we get 



(3.16) 



Since (as before) u>i=i — —y^-u> and k — 8ir, equation 

(|3.15(l becomes the equation (1.3.6), and we thus recover 
the solution (1.3.8). 

Turning now to I > 2 we are confronted with an in- 
triguing problem. Direct calculations (checked and for 
I = 4 done by MATHEMATICA) show that for n = 2 
all functions Q|, Q|, Q\ given by gTJ), (EU and (|3~TT|) 
vanish. In other words, Legendre functions Qy 2 ,, Q3/2' 
Qy* do vanish on the cut! Using the recurrence relation 
(see ^(j; p. 171, with P*f — > Q^, and here we write x 
instead of /1 = cos X ) 



-m+2 



(x) + 2(/i + l) 



(1 



+ (i/-/i)(i/ + p + l)Q£(a;)=0, (3.17) 



we can prove that in fact all 

Qy2 ( x ) = f° r 



I > 2. 



(3.18) 



We are not aware of a mathematical reference where this 
is noticed. 

How then to find regular solutions of the Legendre 
equation I3.1fl for v = |, n = 2 and I > 2? We start 
from the following observations: if solves Legendre 
equation (|3.1|) but for some special values v = vq this 
Q% vanishes, then the function (<9Q^ / 'dv) v=Vfs solves the 
equation for v — v§ and may be non-vanishing. The 
proof can be done easily by a direct calculation. There- 
fore, we wish to calculate dQ l n /dn from the expressions 
(|3~7|) . lO, i|3~TT|) . and put n = 2. Denoting 



<9n 



<9Q 



H4 



f>7? 



(3.19) 



we arrive at the following expressions: 



(I 



1 



sin 5/2 x 



1 



sin(4x) + 2 sin(2x) - 3x , 
(3.20) 



7 



Ql = 



(I 



■ 7/2 

sin ' x 



(3.21) 



sin(5x) 



15 



sin(3x) — 10 sin x + 15% cos x 



2? 



7T\ 2 

2 



sin 9 / 2 X 



- g sin(6x) + 3 sin(4x) 



375 

— sin(2x) - 15x[3cos(2x) + 4] 



(3.22) 



MATHEMATICA yields nice expansions at x -> 0: 



7T\ 2 

2 



/ 7T\ 2 f 



Q 2 

*"(3 



5 21 

|x 7/2 -|x 11/2 + o(x 15/2 ) 



-|x 9/2 -^x 13/2 + o(x 17/2 ) 



(3.23) 
(3.24) 
(3.25) 



Hence, functions H3.2()|I - I|3.22J) are well behaved near the 
origin, 



(const. )x 



i+h 



as x -> 0, 



(3.26) 



however, they at first sight diverge at x = 7r - 

Realizing that x — t is an alternative origin, we take 
as independent solutions (notice that Legendre equation 
is invariant under /i — > —fx) 



(3.27) 



They read 



sin 5/2 x 



1 (4 X )-2sin(2x)-3(7r-x) 



V2 V 1 Si n 7 /2 X 
10 sinx — 15(7r — x) cos X 



I sin(5x) - ^ sin(3x) 



(3.28) 



(3.29) 



375 



sin 9/2 x 



■ sin(6x) — 3 sin(4x) 



• sin(2x) - 15(tt - x )[3 cos(2x) + 4] 
o 

At x — > 7T they behave as Q3.23J) - I|3.25|) with x 
7>2 — (const. )(tt - x)' + ^ as x — » t 1 "- 



(3.30) 

7T, SO 

(3.31) 



The expansions (|3.26l) and l|3.31|l were obtained by di- 
rect calculations for I — 2,3,4. Let us briefly indicate a 
more general argument demonstrating their form for gen- 
eral I. In ^3] one finds (see p. 178) the recursion relations 
(here again x = cos x) 



(3.32) 



and 



(1 - ^ 2 )^ = iy + l)xQ$ - (y - m + l)Qt +1 . (3.33) 

Expressing from the last equation, substituting into 
()3.32(l . taking <9/<9^ and putting v = = I + | and 
a; = cos x, one finds the recursion relation 



a+ 2 )c ° SX ^^ 



smx 



d aQ^ +l 



3/2 



dx dv 



(Z + 2)(Z-2)sin X 



dQl 



3/2 



3/2 



Let us assume that for small x one can write 



Z+i 

9QI7 2 



dv 



= aiX P , 



3/2 



and that 



dQl 



dv 



oti-\X 



(3.34) 



(3.35) 



(3.36) 



3/2 



is known. Then equation (|3.34|l implies p = I + h and 

(3.37) 



(Z + 2)(/-2) 
m = 1 ai-i. 



l + i 

Continuing to / = 2 and regarding l|3.20|l for Q|, we 
obtain 



Oil 



n\i (l + 2)\{l-2)\ 
21 (2i + 1)1! : 



I > 2. 



(3.38) 



Therefore, the asymptotic expansion l|3.26|l can be writ- 
ten for general I > 2 

Q l 2 = a lX l+ K (3.39) 

where a\ is given by l|3.38ll . It is easy to check that 
for I = 2,3,4 one gets the first terms in the expressions 



Next we need to calculate the Wronskians. Defining 
the Wronskians with respect to x (rather than with re- 
spect to /i = cosx), 



dx dx 



(3.40) 



and substituting from lj3~2l1|l - lj3~2^|) and $£^ - $£5$ , 

MATHEMATICA yields: 



W 2 



12tt- 
sinx 



60tt 2 
sinx' 



W 4 



720tt 2 
sinx 



(3.41) 



Of course, the simple dependence on x could have been 
obtained easily since the Wronskians satisfy a simple dif- 
ferential equation. MATHEMATICA is useful to deter- 
mine the coefficients. (If the Wronskians are defined in 
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terms of d/dcosx, the last expressions are multiplied by 
factor — 1/sinx-) Playing with recurrent relations given 
in [13 (p. 171) for both Q£(je) and Q£(-x), and taking 
appropriate derivatives d/dv, we can find the general for- 



l-h, 



mula for Wi, I > 2, in terms of Q 2 ,V 2 and Q 3 j 2 (—x 



Wi= - (Z + 2)(/-2)(e 



Q^ 1 



V l 2 Q 3/ 2 (x) + Q l 2 Q 3/ 2 (-x) . (3.42) 



Since Q 3 / 2 = for Z > 2, the second term in square 
brackets does not vanish only for I = 2. 

As before, we can now solve the inhomogeneous equa- 
tion (|1 . 1 f>|> with k = +1, I = 2,3, ••• by variation of 
the parameters. Remembering that Q l 2 is well-behaved 
at x — 0, whereas P| at X — we write the solution in 
the form 




&2 



(ST%d X ' 



& 



o Wi sin 2 x' 



Wi sin 2 x' 



(3.43) 



where for I = 2, 3, 4, functions Q's and P's are given 
by E2nj-CI22I> and fH3HT3Hll . and (W; sin^ xT 1 = 
(const.) sin 5 x' as a consequence of (I3.41|l . For a general 

14- — 14- — 

I, Q\ = (9Qu 2 /dv) v =3/2, Qv 2 being determined by 
andV l 2 (x) = Q 1 2 (tt-x)- The final physical solution 
is then obtained from w; = (sinx)~ 3 ^ 2 w;, with Wi given 
by H3.43fl • It is easy to make sure that for £1 — ui regular 
at the poles (i.e., due to (5T°)i ~ sin 2 x at the poles), uji 
is well-behaved at the poles. 

Considering (6T®)i concentrated around the equator 
at x € (tt/2 ± A), for example, then near the pole x = 
we find (I > 2) 



uji = 2naix 



i-i 



'-^(ST^dx', (3.44) 

-A Wi sin 2 x 



X we find 



S ™ = "(i) i^-^ + Oix 7 ^), (3.45) 

Q ™ = ^D'f^ + ^ + ^X 772 ), (3-46) 

S " = {iy^-^ 2 + *)x 5/2 + 0(x 9/2 ), (3.47) 

Q 2 W = i(f)^(^ 4 + 5iV 2 + 4) X 5 / 2 +0( X 9 / 2 ),(3.48) 

Ql = - (|) l T ^(-n 6 + 14n 4 -49n 2 + 36)x 7/2 + 
+0(X U/2 ), (3-49) 
QiN = -z(|)^(7V 6 + 14iV 4 + 49iV 2 + 36)x 7/2 + 

(3.50) 



Qn = - f^'-^-(n 8 -30n 6 + 273n 4 - 820n 2 
" V 2 / 945 



+576) X 9 / 2 + 0(x 13/2 ), 



(3.51) 



Q-n = -i 



1 AT 

f EV (jV 8 + 30AT 6 + 273JV 4 + 

\2/ 945 



820i\r 



+576)x 9/2 + 0(x 13/2 )- 



(3.52) 



To obtain these expansions we could not have used for- 
mulas given, for example, in [T(il | (see p. 196), since the 
leading term in Q^(x) at x = 1 is proportional to cos(7r/i); 
hence, it vanishes for /i = 1+ ^ . However, we may proceed 
similarly to the way how we obtained the coefficients ai in 
(|3.38l) . We substitute for Q^ +1 (x) in (|3.32l) from equation 

(I3.33|) . put x = cosx, fx = l + h, ^ = n— i = Q^, 

Tl 2 

and so arrive at the recursion relation 
1 dQi 



Ql— ^+tan X 



(l + n) tan xQ^ 1 = 



' n — I 



n — I dx 



(3.53) 

Regarding results (|3.45|) - (|3.52ll . assume that for small x 



Q n = ai,nX P , 



(3.54) 



and that Q l n 1 = ai-i yU x l 2 is known. Then (|3.53|) im- 
plies p = I + 7j and 



where factor ai is given by 



(n + l)(n- I) 
m ' n = 2TTl ai - hn - 



(3.55) 



Continuing to I — 1 we get 



B. The fluid angular velocity as a source 

To solve 111. 16| l with k = +1 and f2; given, we use 
solutions 13.211 for general v = n — \. In particular, for 
I = 1, • • • , 4 we have explicit expressions H3.5fl - I|3.12|) f° r 
both cases of n real and purely imaginary n — iN. These 
solutions are well-behaved at the origin. Indeed at small 



(-l)'(o) 



For n = iN, 



(XL n = (3.56) 
7T. 1 n(n 2 - l 2 )[n 2 -{I- l) 2 ] ■ • ■ (n 2 - 4)(n 2 - 1) 



(2Z + 1)!! 



Q\ N = ai,iNX l+ K 



(3.57) 
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where 



and 



OLIAN 



N(N 2 + l 2 )[N 2 + (/-l) 2 ] 



(A 2 



(3.58) 
-4)(A 2 + 1) 



(2Z + 1)!! 

The second (independent) solutions which are well- 
behaved at the other pole, x = are, as before, obtained 
by taking Q's at ir — x- Thus, 

K(x) = Q l n(K-X) , -P l iN (x) = Q l iN (n- X ) (3.59) 

are solutions to the homogeneous Legendre equation (|H . 1|) 
on the cut which are regular at x — 7r - 

We do not write them down explicitly since they follow 
easily from Q's given in (|3.5fl - (|3.12f) . At x —> n 



iN 



What are the Wronskians? We get 

sin(n7r) 



m, N {Qj N ,v 2 N } 



n(n — 1)- 



A(A 2 + 1) 



sinx 

sinh(7V?r) 
sinx ' 
sin(n7r) 



n(n 4 - 5n 2 + 4) 
2 smx 



|a(A 4 



5A 2 
sinh(A7r) 



4) x 



smx 



(3.60) 

(3.61) 
(3.62) 
(3.63) 

(3.64) 



where we restricted ourselves to I — 1,2 since the poly- 
nomials depending on n (resp. A) are proportional to 
those in the expansions (|3.45|I - H3.52|I at x ~~ * for 
which we succeeded to derive the formulas (|3.57f) - (|3.59f) 
for general I. (This is due to the fact that the functional 
dependence of the Wronskian is governed by a simple 
differential equation, and the proportionality factor can 
be determined at any point, e.g., at x - * 0.) The gen- 
eral form of the Wronskians can be derived by expressing 
functions V l n in terms of Q^'s and P^'s. Using the rela- 
tion giving Q%(— x) in terms of Q^(x) and Pjfix) [not to 
be confused with T>i = Q£(n - x)] ~ see e.g. 0j, p. 170 



and putting fj, = I + -1 , v = 



we get 



K = q;:a(-*) = (-i)' +1 x 



(3.65) 

Q l n-l ( X ) COS ( n7T ) + \ pl n-\ ( 2; ) Sin ( n7r ! 

Employing then the well known formula (|10|. p. 170) 



W[P£,Q!l] = 
we arrive at 



i 



r(i + z/-/z) i 



(3.66) 



x[n z -{l-l) 2 ] •■■ (n 2 -4)(n 2 -l) 



l+i^n(n a - f 2 ) x (3.67) 
sin(n7r) 



smx 



WlAQmi'PiN] = ^N(N 2 + l 2 )x 



x[N 2 + if 



(A 2 + 4) (TV 2 + 1) x 



(3.68) 

sinh(TVTr) 
sinx 



It is easy to regain <|3.61fHl|3.64f) for I = 1, 2. 

Now we can finally solve the inhomogeneous equation 
p. 16(1 . with ill as a source, again by the variation of the 
parameters. Regarding the relation (|1.14|) between ZJi 
and coi we arrive at the solutions in the form 



n 2 -A |„ ; T (sinx') 3/2 Q , w 



(sinx) 3 / 2 \ n J W t 
(sinx') 3/2 



Q 1 



V^dx' , 



(3.69) 



analogously for n = iN: 



N 2 



(sinx) 3 / 2 



+ Q\n 



V. 



N 



(sinx') 



A3/2 



l,N 



m. 



N 



(3.70) 



All quantities entering l|3.69fl and (|3.70|) are given above 
(for I = 1, • • • ,4 in the explicit forms), the time- 
dependent 'constants' n, N are determined by l|1.17ll . 
which was used to express A 2 in terms of n 2 and A 2 . 
Notice that solutions (|3.70|) are all real since Q\ N and 
V\ N are purely imaginary. We could have divided all so- 
lutions of the homogeneous equation given in the present 
section by n (resp. A) as we did for I = 1 in Paper I, and 
so define *Q l n = Q} n /n, *Q\ N = Q\ N /N etc which, as can 
easily be seen, remain non-vanishing for n ~ 0, N = 0. 
However, the Wronskians l|3.67|l and l|3.68|l are propor- 
tional to n 2 (resp. A 2 ) as n — > (resp. A — > 0), so the 
expressions Q^T^/Wi,™ (the same with A) remain non- 
vanishing in this limit. Hence physical solutions l|3.69fl . 
(|3.70l) for u>i do not vanish when n = (resp. A = 0). 

Now the final solutions l|3.69|l . (|3.70|) for any I can be 
applied to the special situations considered for I = 1 in 
Paper I. So, for large A (I3.68f) implies 



r2i+l 



,JVti 



H"/..y = —A . . 
2 2 smx 



(3.71) 



and the inspection of H3.5fl - (|3.12f) and the relation (|3.60|) 
lead to 



T>\n = *'(-!)' 



7T\f 1 N l e N ^-^ 



2 I 2 sin i/2 x 



(3.72) 



For A large but x small expansions 13.45fl - H3.52fl and 
general formulas l|3.57fl . (|3.59fl give 



(2/ + 1)!! 



(3.73) 
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Hence our solution l|3.70[l for N large near the origin is takes the form 



A'+V" 1 
(2/ + 1)1! 



^(xOsin^'e-^'dx'. (3.74) 



For I = 1 we recover the leading term (independent of x) 
in the expression (1.4.28). The next term in front of the 
integral in (|3.74l) . which is ~ x l+1 > comes from the second 
terms in the expansions (|3.45() - (|3.52l) indicated only by 
C(x i+ 5); their exact forms, if needed, can be obtained 
without difficulty by MATHEMATICA. The resulting w ; 
at given small x can become large at high l's since A 2 = 
Ka 2 (p + p) — 4 may become very large in, say, radiation 



dominated universes (when (p + p) 



) close to big 



bang. Near the perturbation with f2; concentrated in 
r (xo) ± A, with A A < 1, N large, the solution (|3~?0"|) 
leads to 



w;(xo) 



Ne 



-n\ x '- 



smx 



dx', 



,sm xo, 

,(3.75) 

i.e., the result identical for all Z, and thus coinciding for 
I = 1 with (1.4.29), in which A 2 = N 2 + 4, and only the 
leading term in N is kept. 



IV. AXISYMMETRIC TOROIDAL 
PERTURBATIONS FOR k = -1 



With k = — 1, the homogeneous differential equation 
corresponding to the equation l|1.16|) with f2; as a source 



d_ 

Of i 



(1 



v(y + l)- 



W X = 0, (4.1) 



again the Legendre equation. Now \i = \J\ + r 2 = 
coshx, t = sinhx, so the independent variable fi € 
(l,+oo). The solutions of the Legendre equation for /i 
such that its real part is > 1 (fx is then often written as z) 
are sometimes denoted by capital Gothic script letters. 
This is the case in 10] from which we shall again take 
some formulas but we keep the Latin script. If v = |, 
equation (|4.1|) corresponds to the inhomogeneous equa- 
tion (|1.15|) with (ST®) i as a source. 



As in the case k — +1, we start from the general ex- 
pression for the solution in terms of the combination of 
two hypergeometric functions. We first consider the Leg- 
endre function P{f(z) as given in [13 (p. 154, the fourth 



coshx, put fx = I + 



line), substitute for z = fi 
and, as before, v + | = n. Regarding standard prop- 
erties of the hypergeometric functions for special values 
of parameters, we find after simple algebra that func- 

tions P _\ (coshx) can be rewritten as the following fi- 
nite sums: 



p i+i = i r(n + / + i) i Jv M-0fc(* + i)fc 

*-* r(n+l) (sinhx)5\^ {n+l) k k\ 

(cf. |j!OJ0. We shall also use 

p-dH) = i_r(n-Q i J <(-i) k (i + i) k 



/2n T(n + 1) (sinhx) 5 



k=0 



(n + l) fe fc! 



(2 sinhx)' 



1 

(2 sinhx)' 



(4.2) 



(4.3) 



which is also the solution since the Legendre equation use d to find a suitable form of the solutions Q l ^K ■ 
is invariant under the change I + | — > — (/ + |). This 2 

solution can be written as a linear combination of P \ 



and of another solution, Q n _ 2 i, as follows |12fl 



z+i . fK\h (n + l)(n + l- l)---(n+l) 

*n-l - U2 



Q„_A = Mn) " " " 7TT h 

(sum x) ; 



y^ (-l)'+ fc H) fc (/ + l) fc e-( fc +")x 
^ (n + l)&fc! (2sinhx) fe ' ' 



fc=0 



Another useful expression for P_^_l is provided by the 
The last equation combined with (|4.2|) and (|4.3|l can be inversion of Whipple's formula (see e.g. 01 > P- 164) 
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which in our case of \x and v reads 



P 



Q 2 - l)ie m *Q- n (z) 
^~/2 T{1 -n + 1) 



(4.6) 



where z = coshx (notice that the orders and degrees in 
P and Q are here interchanged). 

The case k = — 1 is simpler than that with k = +1. 
The "constant" parameter n = v + 1/2 does not become 
imaginary (cf. equation l|1.17|) l. and the particular case 
v = 3/2, n = 2, corresponding to the homogeneous equa- 
tion in i|1.15|l with ST® as a source, does not require as 
special procedure as for k = +1. The relevant indepen- 
dent solutions of the homogeneous equations are 



V 



P 



4+\ 



(4.7) 



where we use the notation similar to that with k = +1 
but put bars above P's and Q's to distinguish the two 
cases. The solutions (|4.7ll are given as finite sums by 
jO} , respectively 

The asymptotic behavior of Q n and V n at the ori- 
gin and at infinity can be read off from literature where 



asymptotic forms of Pg and Q£ are given (e.g. 01 > P- 
196 and 197). After arrangements we find at small x 



x. 



at large x 



V 

• r 



(-D-d 



(2coshx) r 



(2Z + 1)!! ' 
fir^ (2Z - 1)1! 



7r (n + l){n + 1 - 1) •••n' 
(n + Z)0 + Z-l)' 



(-1) 



(n + 1) 



(2 coshx) 



n+i 



(4.8) 
(4.9) 



(4.10) 
(4.11) 



Before we consider the solutions of the inhomogeneous 

equations we list explicit expressions for V n and Q n for 
I = 1,2,3 as they follow, after suitable arrangements, 
from general formulas l|4.3|l . (|4.5(l and the definition Q4.7[l . 
They look as follows: 



v 1 



fir. 



Or. 



1 



7T / 2n(n 2 — l)(sinhx) 3 / 
1 



- j(n - 1) sinh[(n + l)x] — (n + 1) sinh[(n - l)x]|, 



2(sinhx) 3/2 



V = 



1 



{n + \)e- {n - 1)x ~ (n - l)e- {n+1)x , 

(n + 2) sinhx[— 3 cosh{(n + l)x} 



ir J n(n 2 — l)(n 2 — 4)(sinhx) 5 / 2 

+ (n + 1) hinh \ sinhf /; \ ! + ■'■> sinb[( /< + 2)\ }>. 

1 



2(sinhx) 5/2 



(n + 2) sinhx (n + 1) sinhxe"" x + 3e" (n+1)x 



3e 



-("+2)X 



2\ 3 



(ri + 3) sinhx 15 cosh{(n + 2)x} 



7r / n(n 2 — l)(n 2 — 4)(n 2 — 9) (sinhx) 7 / 2 
(n + 2) sinhx [(n + 1) cosh(nx) sinh(nx) — 6sinh{(n + l)x}] — 15 sinh{(n + 3)x} 



(I) 2 (^xF^l (n + 3)sinhx 



15e -(»+2)x + (n + 2) sinhx X 



x [(n + l)sinhxe" nx + 6e- (n+1)x 



+ 15e-( n+3 ) x 



(4.12) 
(4.13) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 
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A. The angular momentum as a source of ui 

We need to solve the homogeneous equation 1)4.1(1 for 
v = |, n = 2. At first sight we again encounter a peculiar 
case as with k — +1, since T(n — I) in ((4.3|) has a pole 
at I = 2 for n = 2, the expression (|4.4|) appears to have 
also a pole in this case, and explicit formulas (|4.14|) and 
(|4.16(l seem to diverge for n — 2. However, Whipple's 

formula 1(4.6(1 for V n and the following formula 



V l n {z) 



1 



1 



/^72r(; + n+l)(5 2 -l) 1 /4 

_(4.18) 

which, after employing the relation between Q ; " and 



Qi llOj - is its consequence, show no peculiar behavior 
for n = 2. In fact, the last formula can be used to ob- 
tain simple explicit expressions for V n because (e.g. 0, 
formula 8.6.7) 



Qf(z) = (z 2 - 1) 



d 2 Qi 
dz 2 



(4.19) 



where 



n 



i 



i 



[l20x coshx 



(sinhx) 7/2 480' 
-80 sinhx - 15 sinh(3x) + sinh(5x)] , (4.26) 

' f - • ■ ' " (4.27) 



(?)" 



:{5sinhx[e 



-4x 



(sinh x) 7 / 2 

+4sinhx (3sinhxe^ 2x + 6e~ 3x ) ] + 15e~ 5 * j. 

Starting from the well-known formula (e.g. [Tof . p. 165) 
for the Wronskian W{P-",Q£} = e ra ^/(l - z 2 ), z = 
coshx, which is independent of v and, in our case with 
fx = I + |, the value of / changes only the sign of W, we 
easily get 



m,n{Pn, QJ = W*{P 2 , Q 2 } = 



(4.28) 



Employing again the method of variation of parame- 
ters and regarding the asymptotic behavior of V n and 



Q n , we arrive at the well-behaved solutions of the inho- 
mogeneous equation 1(1.16(1 - and hence of the original 
equation i|1.14[> - in the form 



ip ; (z)log^- 



1 



(4.20) 



Wi-i are polynomials of degree I — 1 given in terms of 
Legendre polynomials Pi(z). 

The same result can be derived directly from the origi- 
nal expressions ((4.31) and (|4.5|) by taking the limit n — > 2. 
So, we take as our solutions of the homogeneous equa- 
tions l|4.1fl for u = 3/2, n = 2 the functions 



lim V r 

n^2 



q 2 



lim Q, 



(4.21) 



where P^ and are given by (|4.T|> . (|4.3|) and (|4.B|I . The 

behavior of P 2 and Q 2 at X ~~ * 00 is given by 14.10(1 and 

((4.11|) with n = 2; the first term in the expansions of V n 

and Q l n at x ~^ 0, equations 14.8(1 and 1(4.9(1 . does not 

depend on n, so P 2 at x — * behaves as V n . 

The explicit expressions for I — 1,2,3 follow easily from 
taking the limits n — > 2 in 1(4.12(1 - 1(4.17(1 . One finds 



s 2 



Pi 



Q 2 



^ 2 i(sinh X ) 3 / 2 , 

IT I 6 



(4.22) 



1 



(sinh x) 3// 



2 [coshx(l - 2 sinh 2 x) 



+2 sinh 3 x] , 
'2\ 3 1 



1 



7r y (sinhx) 5 / 2 96 
sinh(4x)] , 

7T\ I 3 

2/ (sinhx) 5 / 2 ' 



[I2x- 8sinh(2x) 



(4.23) 

) + 
(4.24) 

(4.25) 



2 k 



(sinhx) 3 ^ 2 



q' 2 



Wi (sinhx') 2 



x Wi (sinh X ') 2 



(4.29) 



For I = 1 the angular momentum in (0, X ) is given 
by equation ((3.16|) . with sinx' — > sinhx'; w i and {8T®)\ 
are connected with u> and <5T° as in the cases k = 0, +1. 

Noticing that W defined in (1.3.13), is related to Q 2 given 

9^1 



in (|4"2"3|) by W = - y/^j smh x ) ~ 3/2 Q 2 , we easily find 
that for I = 1 the solution 14.2911 goes over into the solu- 
tion (1.3.14). 

Next consider {5T®)i concentrated at x G (xo ± A), 
Xo > 0. Then near the origin x - ► 



.z-i 



2nai\ 
where (see (14.811 



Xo+A 



Xo-A 



Oil 



Ca(-l) 



^(sinhx'^Odx', 
(4.30) 



1 



(2/ + 1)1! 



(4.31) 



B. The fluid angular velocity as a source 

To solve the inhomogeneous equation 1(1.16(1 with k = 

— 1 and Q,i given, we use solutions V n and Q n determined 
by l|4*7|) . (|Oj) and (|4~5|) . For I = 1, 2 ,3 th ey are given 
in explicit, arranged forms by ((4. 12(1 - 1(4.17(1 . The solu- 
tion uji, obtained by variation of the parameters, leads, 



13 



regarding l|1.15fl . to the following physical solution well- 
behaved at the origin and vanishing at infinity: 



(sinhx) 3 / 2 



-A 2 n t )(sinhx / ) 3/2 



_^ [°° QL(-A 2 ^)(sinh X Q 3 / 2 dy 



(4.32) 



Substituting for W;,„ from (|4^S|) and for X 2 (t) = n 2 (t) - 
4 as it follows from H1.17fl with k = — 1, we finally have 



'L^°°QLa(sinh X ') 5/2 dxj. 



(4.33) 



For I = 1 we notice that with and 7^ we regain the 
solution (1.4.34). 

Assume J7; non-vanishing at large \- Then the solu- 
tion H4.33fl and the asymptotic results for Q n and V n , 
equations Q4.8|) - Q4.11[) . imply that at small x 

(n 2 -4)(n + l)(n + l-l) ■ ■ ■ (n + 1) 

U\ = " T7Z~, ; 7771 " "X x 



4(2? + 1)!! 
x / (sinhx') 2 e _ " x '^dx'. 



(4.34) 



Close to the perturbation, 



LU, 



1 {n 2 

2 ? 



4) 



/ sinh 
o \ sinn Xo 



e -n|x'-Xo| fi; ^ 



(4.35) 

Here I does enter the result directly. So for 2 = 1 it goes 
over immediately to (1.4.38). It is also easy to see that 
(I4.34|) gives for I = 1 the first terms in (1.4.37). 
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